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For fifty years, cosmic-ray air showers have been detected by their radio emission. We present 
the first laboratory measurements that validate electrodynamics simulations used in air shower 
modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from 
charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air 
shower. This experiment provides a suite of controlled laboratory measurements to compare to 
particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air 
shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular 
distribution, polarization, and spectral content. In particular, we confirm modern predictions that 
the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle 
and show that the simulations reproduce the data within systematic uncertainties. 

PACS numbers: 95.55.Vj, 98.70.Sa, 29.27.-a 

The highest energy cosmic rays arrive at Earth with 
energies in excess of 10 2 ° eV. Despite decades of work 
meant to uncover their sources, their origin remains elu¬ 
sive. Observations are limited by the low flux at the end 
of the cosmic-ray spectrum, calling for the development 
of new techniques with high duty cycles, high precision, 
and large surface areas. One promising technique makes 
a measurement of the radio-frequency electric field from 
a cosmic-ray air shower, which is nearly linear with the 
energy of the primary particle. 

Radio emission arises from a cascade of charges moving 
inside a dielectric and in the presence of magnetic field in 
two main ways. Askaryan radiation forms from a charge 
excess built up in the shower due to Compton, Bhabha, 
and M0ller scattering and positron absorption, forming a 
current along the shower axis and radio frequency emis¬ 
sion p] . Such emission has been measured in accelerator 
experiments nu. Geomagnetic emission forms when the 
Lorentz force acts on charges in the shower, generating 


a time-varying transverse current. The former is present 
even with no magnetic field, and the latter is present 
even without a charge asymmetry. In practice, experi¬ 
ments detect the sum of these two effects. 

Several experiments have detected radio emission from 
cosmic-ray air showers [ZH22]- To progress from event 
detection to measurement of a differential energy flux, 
one needs to know the intensity and angular distribu¬ 
tion of the radio emission and its frequency dependence. 
To date, predictions of these parameters have relied on 
simulations and measurements of air showers themselves, 
which can be subject to uncertainties in geometry and 
hadronic interactions in addition to uncertainties related 
to the radio emission. In contrast, this work provides a 
direct laboratory benchmark for the simulations of the 
radio emission using particle showers with well-known 
cascade physics that develop in a precisely known target 
and geometry. 

Our particle showers developed in a dense plastic in- 
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FIG. 1. (a) Schematic of the experiment, not to scale, (b) Simulated (Endpoints) electric field, \E\, in the x-z plane with full 
magnetic field, 131 pC, and 4.35 GeV. 


stead of rarefied atmosphere. Increasing the density, p, 
decreases the shower length, thereby reducing the total 
amplitude of both radiations integrated over the whole 
shower. The total Askaryan radiation scales with the 
shower length as 1/p, weighted by the charge asymmetry. 
Because the magnetic radiation from a given position in 
the shower is proportional to the transverse drift veloc¬ 
ity of electrons and positrons, the magnetic radiation per 
unit track length scales with the magnetic field B and the 
scattering length [23] . or B/p. The total magnetic radi¬ 
ation for the entire shower scales as 5/p 2 , meaning that 
the expected ratio of magnetic to Askaryan radiations 
scales as B/p. So, in order to achieve the same relative 
emission amplitude found in air showers, we increased 
the magnetic field from the terrestrial 0.5 Gauss to 1000 
Gauss, commensurate with an increase in density from 
air to plastic. 

Two formalisms, ZHS m and Endpoints m , are 
at the heart of recent simulations of radio emission, 
ZHAireS [26. and CoREAS [27], respectively. They both 
treat each shower particle track as an independent ra¬ 
diator, summing up the emission from all tracks in the 
cascade to obtain the signal that would be received by an 
observer. The ZHS technique has been adapted to the 
time-domain [28] and calculates the vector potential of 
each particle track. The Endpoints formalism sums the 
radiation due to the acceleration between discrete seg¬ 
ments of the particles’ trajectories [25]. Both formalisms 
are compared to the accelerator data in this paper. 

The T-510 experiment, shown schematically in Fig.[l£i, 
took place at the End Station Test Beam (ESTB) in End 
Station A (ESA) at SLAC. Bunches of electrons with 
energy 4.35 or 4.55 GeV passed through 2.3 radiation 
lengths of lead pre-shower and entered a plastic target, 
generating compact showers with total energy equivalent 
to a ^ 4 x 10 18 eV primary cosmic ray. 

The target was formed from 5.08 cm x 10.16 cm x 
30.48 cm bricks of high-density polyethylene (HDPE) 


with density 0.97 g/cm 3 . Being 4 m long, 0.96 m tall, 
and 0.60 m wide, it contained the majority of the parti¬ 
cle shower. The beam was 0.13 m above the bottom of 
the target. The bricks on the top surface were machined 
to a 9.8 degree angle below horizontal to avoid total inter¬ 
nal reflection. The index of refraction of HDPE is 1.53, 
resulting in a refracted Cherenkov angle of 28° from the 
horizontal (49° before refraction). Rays emanating from 
the target at the refracted Cherenkov angle intersected 
the antenna plane at 6.5 m above the position of the 
beam. The target floor was lined with an RF absorbing 
blanket. Several pieces of ANW-79 absorber were placed 
at both sides of the target and at the exit surface of the 
target. 

Fifteen water-cooled coils placed under the target were 
used to create a vertical magnetic field of up to ±970 G. 
This was achieved by supplying sets of five coils in se¬ 
ries up to 2400 A direct current with reversible polarity. 
The coils were aligned along the beam axis on a 10.16 cm 
thick steel plate and were staggered in two rows in order 
to create a more uniform magnetic field. The vertical 
magnetic field measured at the beam height and maxi¬ 
mum current had a RMS variation of 72 G. The magnetic 
field was strong enough to bring the expected intensity 
from the magnetic effect to the same order of magnitude 
as that expected from the Askaryan effect. 

Four dual-polarization, quad-ridged horn antennas 
used in the ANITA experiment [29] recorded the elec¬ 
tric fields generated in the particle shower. An overhead 
crane allowed for the movement of the antenna array to 
sample the electric field at many positions. The antenna 
array was placed at the far wall of ESA, 13.5 m from the 
entrance of the beam to the target. It was tilted at 19.6° 
to better align the expected radiation with the boresight 
of the antennas. The antenna array covered vertical dis¬ 
tances between 1.4 m and 12.4 m, corresponding to angles 
between 40° and 55° with respect to the beam line at the 
beam entry point, which, due to the lead pre-shower, is 
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close to shower maximum. 

The antennas used are sensitive to the 200-1200 MHz 
band. The comparable frequency range in air showers 
is lower, because it is inversely proportional to both the 
Moliere radius and sin#, where # is the observation an¬ 
gle. The Moilere radius scales as 1/p and sin# scales as 
y/p to first order. Taken together, we expect that the fre¬ 
quencies scale as yfp. The T-510 bandwidth translates to 
approximately 10-60 MHz in air showers, comparable to 
the bandwidth of ground-based air shower experiments. 

Signals from each antenna ran through 15.24 m of 
LMR240 coaxial cable and a low-pass filter with a 3 dB 
point of 1250 MHz to avoid aliasing during data acquisi¬ 
tion. Time-series voltages from the horns were collected 
on 2.0 GHz, 5 GSa/s oscilloscopes. A global trigger was 
provided by the broadband transition radiation produced 
by the beam exiting the beam pipe, collected in an S- 
band horn antenna. Events were recorded at 1 Hz. 

The beam charge was measured using an integrating 
charge transformer situated between the beam pipe and 
the target. The mean bunch charge was 131 pC, with a 
shot-to-shot standard deviation of 3 pC. Measurements 
of the bunch charge at several positions indicate a 2% 
systematic uncertainty. 

Particle showers were simulated with GEANT4, us¬ 
ing the measured magnetic field. Five thousand primary 
electrons were injected, and the results were scaled to 
131 pC. The radio emission was then simulated following 
the ZHS and Endpoints formalisms. In calculating the 
radiation from each track we included refraction, Fres¬ 
nel coefficients, and demagnification effects m at the 
surface of the target. 

The simulation, shown on the right in Fig. [TJd, demon¬ 
strates that the expected radiation forms a ring when 
projected onto a two-dimensional plane 13.5 m from the 
entrance to the target, peaking at about 6.5 m above the 
shower axis. The electric field strength, |£j, map shows 
the superposition of magnetic and Askaryan components. 
Since the magnetic field is vertical at the shower, the for¬ 
mer is horizontally polarized, whereas the polarization 
of the Askaryan contribution points radially from the 
shower axis. The interference between the two produces 
the left-right asymmetry shown in the figure. Refraction 
at the target surface makes the ring elliptical rather than 
circular. The ring is cut off on both sides due to the fi¬ 
nite target. The simulation is done using ray optics. The 
top of the target acts as a diffractive slit, with a Fresnel 
zone of about 60 cm at 300 MHz, which is smaller than 
the length of the target. The target width (±30 cm) cor¬ 
responds to a phase lag of about 30°, and so we expect 
the simulation to modestly overestimate \E\ at low fre¬ 
quency. Reflections were not included in Fig. [TJd, but are 
discussed below. Effects due to transition radiation were 
estimated and found to be two orders of magnitude below 
the Askaryan radiation. 

By design, at our antenna locations, the two types 



FIG. 2. (a) Simulated and measured voltages 6.52 m above 
the beam in the horizontal polarization channel at full mag¬ 
netic field, (b) Corresponding power spectral densities. 

of radiation separate into orthogonal polarizations. The 
magnetic is horizontally polarized, while the Askaryan is 
vertical. However, in real antennas, the vertical signal 
leaks into the horizontally polarized channel at about 
the 25% level in amplitude (about —12 dB in power). 
We eliminate this leakage in the horizontal polarization 
by construction in using the difference between field-up, 
Vb + (£), and field-down, Vb_ (£), data, namely V(t) = 
\{Vb+ (t) — VB-(t)). Each waveform is also scaled by the 
beam bunch charge to 131 pC. When this construction 
is applied to the vertically polarized signal, we see only 
noise, as expected. The resulting single-acquisition wave¬ 
form from the horizontal polarization at maximum mag¬ 
netic field strength is shown in Fig. IK Fig. [2)3 shows 
the power spectral density normalized by the sampling 
interval 200 ps/point for the data and 100 ps/point for 
the simulations. 

For comparison with the data, the simulated electric 
fields were convolved with the measured antenna effec¬ 
tive heights [29] and response due to the cables and filters 
used in the system. The convolution was performed in 
the frequency domain following standard techniques m- 
The predicted values in Fig. [2] from the ZHS and End¬ 
points formalisms agree to within 3% in peak amplitude 
and 7% in integrated power. The dominant features in 
the time-domain waveform arise from antenna response 
and filters, which the two simulations have in common. 
The shape of the simulated waveforms reproduces the 
data well, giving us confidence in the experimental mod¬ 
eling. The absolute scale is discussed below. 

Internal reflections from the bottom of the target in¬ 
terfered with the signal transmitted through the top sur¬ 
face of the target, which is apparent in the modulation 
of the power spectral density in Fig. It- A low ampli¬ 
tude reflection with two internal bounces arrived at the 
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FIG. 3. (a) Horizontally polarized signal normalized by ver¬ 
tical showing the expected linear behavior vs. magnetic field, 
(b) The oscilloscope traces (solid) show the polarity flip. Mod¬ 
els (dashed) are shown for opposite polarities. 


antenna ^ 7 ns after the main pulse, and was responsi¬ 
ble for the A/ ~ 150 MHz frequency beat in Fig. H>- 
More important is the radiation that reflected only once 
off the bottom of the target. Assuming that the absorber 
has a higher index of refraction than the HDPE, that re¬ 
flection arrived at the antenna ^ 1 ns after the direct 
pulse, but with inverted polarity. Adding a time delayed 
equal amplitude reflection to the direct pulse for a par¬ 
ticular antenna varies the simulated peak time-domain 
amplitude by a factor of up to 1.38 for horizontal polar¬ 
ization and 1.43 for vertical polarization. Thus, averaged 
over 300-1200 MHz, the uncertainty in the models due to 
reflections is of order 40% for both polarizations. We ex¬ 
clude the 200-300 MHz band in the data and simulation 
comparisons due to uncertainties related to diffraction 
and in the antenna response at low frequencies [29] . For 
the horizontal polarization shown in Fig. [2^i, the time- 
domain peaks of the data exceed the simulations by 35%, 
commensurate with the systematic uncertainty. 

Fig-! shows that the amplitude of the horizontally 
polarized emission is linearly dependent on the magnetic 
field. The polarity of this induced voltage changes sign 
when the direction of the magnetic field flips direction, in¬ 
dicating that the transverse current flows in the opposite 
direction. The vertically polarized emission is observed 
to be constant with respect to magnetic field strength. 
The difference in slopes between the data and simula¬ 
tion is 20% which given our current systematic uncer¬ 
tainty should be taken as agreement between the two. 
This agreement along with the expectation that the ra¬ 
tio of magnetic emission to Askaryan emission scales with 
B/p confirms that transverse currents generate magnetic 
emission in air showers. 


Several aspects of the radio emission in cosmic-ray air 
showers contribute to the formation of a conical beam 
pattern centered around the Cherenkov angle [32, 33 . In 
this experiment, the angular radiation pattern was mea¬ 
sured by placing the antennas at different vertical posi¬ 
tions. The power profile, which traverses the expected 
peak of the cone, is shown in Fig. [4] for three different 
frequency bands. Each profile is normalized by its total 
power. 

The observed cones are somewhat different than those 
from air showers. In both cases, the expected power spec¬ 
trum observed at the Cherenkov angle peaks at a fre¬ 
quency determined by the transverse size of the shower 
Lx, which is a fraction of a radiation length Xq. The 
Cherenkov cones have widths SO = ccj)/ (n/Ly tan Oc), de¬ 
termined by the angle over which the shower is coherent, 
</>, the frequency of observation, /, the Cherenkov angle, 
Qc, and the shower length, Ly. At the peak power fre¬ 
quency, the width of the Cherenkov cone is determined 
simply by the aspect ratio of the shower SO = Lj_/Ly. 
For the T-510 beam of 4.5 GeV electrons, Ly ~ 2Xo and 
we both predict and observe SO ~ 5°. For air showers 
induced by 10 17 eV primaries, Ly 5Xo, and the in¬ 
ner edge of the Cherenkov cone is washed out because 
Oc (i°) < S0( 2°), causing the Cherenkov feature to ap¬ 
pear as a filled-in disk, as observed by LOPES [341. At 
higher frequencies, the width of the Cherenkov ring scales 
as 1//, and the ring becomes well defined, as observed at 
LOFAR [22: and inferred by ANITA [35] . The accurate 
simulation of the Cherenkov cone and frequency behavior 
shown in Fig. [4] is directly relevant to the ability of the 
simulations to model emission from air showers. 

The signal polarization observed in T-510 confirms the 
paradigm that transverse currents due to the geomag¬ 
netic effect and longitudinal currents due to charge excess 
produce the radiation observed in air showers, consistent 
with measurements by CODALEMA [36] and AERA [2T] . 
The observed Cherenkov cone in this work and at LOFAR 
[22] indicate that the refractive index is an important 
component to accurate modeling of the electrodynamics. 
Recent comparisons between the microscopic calculations 
and LOPES data confirm that first-principle calculations 
accurately predict the absolute scale of the radio emis¬ 
sion [37] . but the measurement is subject to uncertainties 
associated with air shower physics such as composition 
and hadronic interaction models. In using a fixed beam 
geometry and electromagnetic shower composition, this 
experiment confirms the absolute scale of the microscopic 
calculations with different systematic uncertainties. 

We have presented the first laboratory benchmark of 
radio-frequency radiation from electromagnetic cascades 
under the influence of a magnetic field. We compared 
the radio emission produced in a well-defined target ge¬ 
ometry with a well-defined particle shower to predictions 
made by microscopic models, which rely on first princi¬ 
ples of electrodynamics and have no free parameters. The 
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FIG. 4. Beam patterns for three frequency bands in the hor¬ 
izontal polarization at 970 G. 


models agree and accurately predict the absolute scale of 
the radio emission to within our systematic uncertainty. 
The observed radiation grows linearly with magnetic field 
strength and forms a conical beam centered around the 
Cherenkov angle. Being in a relevant frequency range and 
independent of hadronic interaction models and complex 
geometries, this experiment is complementary to in situ 
observations of radio emission from air showers by cur¬ 
rent experiments mmmmm- It also strengthens 
the case for proposed experiments based on the radio 
technique [38MTL T-510 provides strong evidence that 
the electromagnetic simulations can be used to reliably 
predict the radio emission from extensive air showers. 
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FIG. SI: (a) The HDPE target stacked on top of the magnetic held coils. The quad-ridged horn array in the y-z plane (b) 
shown close to the target for scale and (c) in a position for data-taking (12.4 m above the beam, and 13.5 m from the beam 
entrance to the target.) 


Supplemental Material: Accelerator measurements of magnetically-induced radio 
emission from particle cascades with applications to cosmic-ray air showers 


I. THE T-510 EXPERIMENTAL DESIGN 

The T-510 experiment was designed to call out the critical features of geo-magnetic radiation. To do this, we used 
a geometry that separated the magnetic radiation and the Askaryan radiation into two different polarization states. 
Fig. SI and Fig. la show the experimental setup, where an HDPE target was constructed on top of a series of magnetic 
held coils. Radiation produced in the target was recorded by an antenna array shown in the near held in Fig. Sib. 
For the data presented in this work, the antenna array was further away (13.5 m) from the entrance of the beam to 
the target, as shown in Fig. Sic. 

The magnetic held map at maximum current is shown in Fig. S2, measured at the beam height, 21 cm above 
the coils with 5 cm x 5 cm grid spacing and 3.64 G precision. There is a 6% uncertainty on the magnetic held 
measurements arising from the average percent error of the Hall probe used. Because the vertical (B z ) held falls off 
near the edges of the coils, the edge of the target was placed in the middle of the hrst coil. Along the beam position, 
the average magnetic held in the vertical direction is 845 G, and the RMS variation is 72 G. The peak magnetic held 
strength is 970 G. In our simulations of the T-510 experiment, we included the measured, three-dimensional magnetic 
held. 

As shown in Fig. S2, the magnetic held along the beam is strongest in the vertical (z) direction. Because of this, the 
transverse current that develops in the shower is primarily oriented along the transverse (x) axis, thereby producing 
horizontally-polarized radiation. In contrast, the Askaryan radiation that arises from the longitudinal current in the 
shower is radially polarized. When the antenna array is conhned to the line dehned by the beam and the vertical 
direction, the magnetic radiation observed at the antennas is horizontally polarized, while the Askaryan radiation is 
vertically polarized. While the ZHS and Endpoints models calculate the radiation per particle track and are therefore 
ambivalent to macroscopic effects, the simulation shown in Fig. lb confirms this simplified picture. 


II. SCALING TO COSMIC RAY AIR SHOWERS 

The goal of the T-510 was not to accurately replicate all the conditions of an air shower in the lab. Rather, the 
goal was to create a model system with relative field strengths arising from both magnetic and Askaryan emission 
comparable to those found in air showers. Even though direct comparisons between cosmic ray air showers and the 
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FIG. S2: The measured magnetic field at the beam height along the vertical (z), longitudinal (y), and transverse (x) axes. 


showers in T-510 are difficult, the radiations observed in T-510 are relevant for air shower experiments, because the 
densities, frequencies of observation, and the magnetic field strength are scaled appropriately. 

In air showers the geomagnetic emission dominates over Askaryan emission, but in a dense dielectric such as HDPE, 
the Askaryan radiation dominates because it scales with the shower length. The source for both radiations is the 
current distribution in the shower. For Askaryan radiation, the total longitudinal current depends on the total track 
length of particles, weighted by the charge asymmetry. The total track length scales as the radiation length or as 
1/p. For the magnetic radiation, the total transverse current, J± depends on the total track length weighted by the 
transverse velocity, or the accumulated transverse acceleration, a_ j_, during a radiation length. J_\_ ~ a±t over time t. 
The acceleration increases with B , and the time of acceleration increases as the radiation length or 1/p. The total 
transverse current therefore scales as B/p 2 . To recreate a comparable ratio of source currents for magnetic emission 
to Askaryan emission as observed in an air shower, we preserved the ratio (B/p 2 )/(l/p) = B/p. To mimic the relative 
importance of geomagnetic radiation to Askaryan radiation in air showers, the quantity B/p in T-510 is comparable 
to that in air showers. 

Two delay times, the transverse delay time Sti and the longitudinal delay time ^ 2 , determine the frequencies and 
angles over which showers are coherent. The transverse delay time determines the range of frequencies where the 
radiation is in phase, while both the longitudinal and transverse set the range of angles over which the radiation is 
coherent. The bandwidth observed in T-510 is higher than the corresponding bandwidth for air shower experiments. 
The shape of the observed beam pattern in depends on both the Cherenkov angle and beam width, both of which 
differ in T-510 and in air shower experiments. 

We first assume that the radiation is coherent over the transverse dimensions of the shower, L_ l, which we take 
to be one twelfth the Moliere radius, R m . The observation frequency, /, is comparable to the time difference, St\, 
between radiation from the center of the shower (point a) and from the edge of the shower (point b ) at the observation 
angle, #, as shown in Fig. S3. The transverse delay time is then given by 


Sti 


nLj_ sin 0 
c 


(SI) 


where 0 is the observation angle. At the Cherenkov angle, Oc, f ~ 1/St\ ~ c/(Lj_sinOc) ~ c/(R rn sinOc)- In T-510, 
sin 6c = 0.8, but in air showers sin 0c scales as 1/yfp. Therefore, the T-510 measurements made at 200-1200 MHz 
are comparable with observations of air showers in the 10-60 MHz band. 

Assuming that the angle 0 in Fig. S3 is given by 0c + S0, the power spectrum at 0c peaks at a frequency /o = 4>/Sti , 
corresponding to the highest frequency over which the shower is coherent. On the Cherenkov cone, SO = 0, such that: 


h = -r^~r ” t-tV < S2 > 

nL _l sm 0c R±_ tan 0c 

The longitudinal delay time is the delay time across the length of the shower, or from point c relative to point a 
in Fig. S3. It includes the time required for the particle shower to propagate from a to c and the delay time between 
radiation emitted at point a and point c. 
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FIG. S3: Schematic of transverse (Sti) and longitudinal (St 2 ) coherence times. 


St 2 = —-- cos 0 = — (1 — n cos 0) (S3) 

ccc 

In solid media, the width of the Cherenkov cone is small: SO « Oc- 

cos 0 ~ cos 6c — SO sin 6c (S4) 

St 2 — — -SOt&nOc (S5) 

c 

Again, assuming the angle over which the radiation is coherent is (j) = fSt 2 = /SO tan 6c , then the width of the 

Cherenkov cone is: 


SO = C( ^ C( ^ /gg\ 

fnL\\tzm6 c f f 0 nL\\ tan 0 C fL\\ 

At / = /o, the width of the cone is determined by the shape of the shower. The width is narrower at higher 
frequencies and broader at lower frequencies. From the simulations, the showers in T-510 are two radiation lengths 
long and the cutoff frequency is approximately 1 GHz, and so SO -5°, similar to the observed widths in Fig. 7. The 
beams observed in T-510 are narrow cones, because 6c > SO. 

Finally, we note that the density of the atmosphere is changing in air showers, while it is fixed for T-510. This 
is important when the the radiation length becomes comparable to the atmospheric scale height, which does occur 
for highly inclined showers such as those observed by ANITA [SI] and proposed sub-orbital experiments [S2, S3]. 
However, for the less-inclined showers observed by ground arrays, the density at the point of emission rather than the 
density gradient determines the radiation strength. 


[51] P. W. Gorham et al., Astropart. Phys. 32, 10 (2009). 

[52] P. W. Gorham et al., Astropart. Phys. 35, 242 (2011). 

[53] A. Romero-Wolf et al., arXiv: 1302.1263 (2013). 









